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摘 要 : 基于 带 有 叶轮 转子 的 双 吸 叶 轮 的 概念 ， 设 计 了 一 种 新 型 无 轴 双 吸 式 血液 和 泵 。 为 了 研究 血液 和 泵 的 水 力 性 能 和 溶血 性 能 ， 
对 设计 工 况 下 血液 泵 进行 了 三 维 非 定常 数值 模拟 ， 同 时 考虑 作为 人 类 心脏 用 时 的 搏动 流 工 况 ， 以 此 为 边界 条 件 计算 分 析 了 搏动 
流 工 况 下 的 内 部 流 场 的 速度 、 压 力 和 壁面 剪 切 应 力 分 布 以 及 对 泵 运行 稳定 性 的 影响 。 数 值 计算 结果 表明 ， 捕 动 流 工 况 下， 实现 
泵 流量 与 扬程 的 波动 输出 ， 与 非 捕 动 流 相 比 ， 其 压力 稳定 性 较 差 。 该 结果 为 血液 泵 在 人 工 心脏 、 人 工 透析 等 现代 机 械 辅助 系统 
的 应 用 提供 有 益 的 参考 。 
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Internal Flow Analysis of a Shaft-less Double-suction Blood Pump 
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Abstract: Based on the concept of double-suction impeller combined with motor rotor, a novel shaft-less blood pump has 
been designed. In order to investigate the hydraulic performance and the biocompatibility of the pump, the unsteady 
numerical simulation for three dimensional flows in entire pump passage near the design condition is conducted. In the 
calculations, the pulsatile flow of a natural artery is used as the boundary condition at pump inlet so as to observe the flow 
characteristics such as the velocity, pressure and wall shear stress in flow passage, and the effect on the operating stability 
of the pump. Based on the results of numerical simulation, the internal flow will be further discussed. 
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0 HJ 
血液 泵 是 以 血液 作为 输 运 介质 的 一 类 医用 血液 式 血 液 泵 具有 噪音 低 、 机 械 摩擦 损失 小 及 使 用 寿命 
和 泵 ， 可 代替 人 体 器 官 的 部 分 功能 ， 用 以 辅助 血液 循 长 等 优点 。 而 相对 于 轴 流 泵 ， 离 心 式 血液 泵 具有 转 
环 ， 是 现代 医学 治疗 急 慢 性 器 官 功能 衰竭 的 主要 装 速 低 、 溶 血 破坏 小 等 优点 。 因 此 ， 有 必要 进行 动 压 
置 。 目 前 ， 常 用 的 血液 泵 多 用 机 械 轴 承 作为 转轴 支 式 血 液 泵 的 研究 。 2009 年 罗 先 武 等 人 I 设 计 一 种 叶 
撑 ， 机 械 轴 承 的 摩擦 损失 大 、 易 损坏 ， 且 运行 时 产 轮 外 径 为 40mm 的 无 轴 超 小 型 泵 ， 该 和 泵 的 转子 由 去 
生 一 定 的 热量 ， 若 用 于 血液 循环 系统 ， 易 产生 凝 压 支撑 ， 这 种 无 轴 双 吸血 液 泵 不 会 造成 大 量 血细胞 
和 溶血 。 因 此 ， 近 年 来 越 来 越 多 的 学 者 开始 投入 磁 破坏 。 针 对 血液 泵 内 流 场 的 研究 ， 澳 大 利 亚 麦 考 瑞 
悬浮 、 液 体 动 压 为 支撑 的 血液 泵 研究 50。 液体 动 压 
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大 学 的 Qian Yi 教授 较 早 开 始 了 关于 血液 泵 的 研究 
工作 加， 并 对 液体 压力 悬浮 式 离心 泵 内 部 流动 特 
性 及 液体 压力 支撑 性 能 进行 了 数值 仿真 研究 ; 
Burgreen 等 人 对 带 有 永 磁 文 撑 转子 的 离心 式 
HeartMate III 左 心 辅助 系统 在 设计 条 件 下 的 动力 学 
特性 进行 数值 仿真 和 分 析 工 作 中 ,仿真 计算 与 实验 
结果 具有 良好 的 一 致 性 。Tsukamoto 等 人 [开发 了 

种 带 有 沟 模 轴承 的 血液 汞 ， 该 泵 带 有 一 种 软 性 材 
料 支 撑 液 体 动 压轴 承 。 

本 文 针 对 一 种 叶轮 外 径 为 24mm 的 新 型 无 轴 双 
吸 式 血 液 泵 在 非 捕 动 及 搏动 流 工 况 下 的 内 部 流 场 进 
行 了 三 维 全 流 道 非 定 常数 值 模拟 ， 分 析 两 种 流动 工 
况 下 血液 泵 的 内 部 流动 特性 ， 探 讨 速度 、 压 力 及 壁 
面 剪 切 应 力 的 非 定常 变化 规律 ， 为 搏动 流 工 况 下 
液 泵 的 运行 稳定 性 设计 提供 参考 。 


1 血液 泵 的 结构 设计 

液 泵 的 整体 结构 主要 由 两 大 部 分 组 成 : 驱动 部 
件 和 流动 部 件 ， 如 图 1 所 示 。 了 驱动 部 件 包括 由 电 枢 绕 
组 、 硅 钢 片 及 永 磁 体 组 成 的 电机 和 控制 器 ; 流动 部 件 
包括 叶轮 、 泵 碗 及 出 口 管 等 。 
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图 1 血液 泵 的 整体 结构 方案 

Fig. 1 Schematic diagram of the miniature pump. 
叶轮 转子 和 压 水 室 之 间 设 有 液 膜 间隙 ， 当 泵 内 
流体 流动 时 ， 该 液 膜 间隙 可 以 作为 叶轮 转子 和 泵 过 
之 间 的 压力 膜 ， 对 叶轮 转子 起 支撑 和 润滑 作用 。 液 
体 经 过 叶轮 后 ， 叶 轮 出 口 压力 高 于 叶轮 进口 压力 ， 


在 压力 差 作 用 下 ， 液 体 沿 着 间隙 从 叶轮 出 口 流向 叶 
轮 进口 ， 冲 刷 间隙 流 道 ， 防 止 在 流 道中 形成 流动 死 
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区 ,并 带 走 电流 通过 线圈 及 转子 旋转 所 产生 的 热量 。 
该 方案 可 有 效 防止 由 于 高 温 及 血液 长 时 间 停 滞 和 积 
淤 而 导致 的 凝血 现象 ， 而 采用 液体 动 压 悬 浮 支 撑 ， 
可 避免 由 于 机 械 轴承 对 红细胞 破坏 的 溶血 现象 。 


2 数值 计算 方法 
2.1 计算 模型 与 网 格 划 分 

全 流 道 模型 由 进口 管 、 叶 轮 、 压 水 室 、 出 口 管 四 
部 分 组 成 。 以 叶轮 区 域 和 旋转 吸入 段 (进口 管 中 间 部 
位 ) 作为 旋转 坐标 系 ， 而 进口 管 、 压 水 室 、 出 口 管区 
域 作 为 静止 坐标 系 处 理 。 本 文选 用 SST k- o mR 


型 进行 数值 计算 。 
综合 考虑 数值 计算 精度 和 经 济 性 , 选用 节点 数量 


为 86 万 的 结构 化 网 格 。 ae S ie 
为 了 更 加 精确 地 计算 叶片 表面 的 压力 及 剪 切 

a oa ii an a E 

小 网 格 间距 为 30um。 
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图 2 计算 模型 网 格 划 分 示意 图 
Fig. 2 Mesh generation 
2.2 搏动 流 工 况 及 模拟 方法 
根据 人 体 搏动 流 流 量变 化 曲线 , 取 心 脏 搏 动 周 期 
T=0.5s 。 泵 进口 给 定 流 量变 化 规律 ， 图 3 为 两 个 周 
期 内 泵 的 进口 MEO) 质量 流量 随时 间 的 变化 规 
律 。 
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图 3 泵 入 口 质量 流量 随时 间 变 化 曲线 
Fig. 3 Variation of mass flow at pump inlet 
一 个 搏动 周期 内 质量 流量 随时 间 的 变化 关系 | 
线 可 用 如 下 函数 表示 : 


lin 


Q(t) =a, + > (a, cos(not) + b, sin(nat)) 


n=1 
(n=1,2...) (1) 
式 中 各 参数 值 见 表 1. 


表 1 质量 流量 变化 函数 参数 表 


Table 1: Coefficients in Equation (1). 


ao ay az a3 ay as ag 


0.0160 0.0027 -0.018 -0.000 0.0018 0.0025 0.001 
82 48 54 26 54 73 17 
续 表 1 质量 流量 变化 函数 参数 表 
bi b bs ba bs be w 
0.0190 0.0026 -0.009 -0.002 -0.000 0.0009 12.5 
64 14 36 11 76 36 7 

3 计算 结果 与 分 析 


3.1 非 搏动 流 工 况 下 血液 泵 的 非 定常 流动 模拟 

在 设计 转速 (5000r/min〉 下 ， 对 设计 流量 工 况 
(CO=O ) 下 的 血液 泵 内 的 流动 进行 非 定常 数值 模 
拟 , 分 析 叶 轮 和 压 水 室内 部 的 速度 、 压 力 及 壁面 前 切 
应 力 分 布 。 


定义 如 下 无 量 纲 的 静 压 系数 C, 和 全 压力 系数 


pt? 


P—Po 


C.= (2) 
? 0.5pU; 

P- Po 

= (3) 
P" 0.5 U2 


式 中 p HHE (Pa) ， pi 为 全 压 Pa), pN 
流体 密度 ， py 为 参考 压力 ， 这 里 取 叶 轮 进口 的 平均 
静 压 


为 了 观察 叶片 旋转 对 内 部 流 态 的 影响 , 将 一 个 周 
期 分 成 5 个 叶片 位 置 对 应 的 时 刻 , 图 中 黑色 轮廓 线 表 
示 同 一 个 叶片 .t=0 与 1=7T 时 刻 对 应 叶片 起 始 位 置 ， 
及 旋转 一 周 后 的 位 置 。 为 了 便于 分 析 叶 轮 及 压 水 室内 
部 的 流动 情况 ， 对 叶轮 叶片 及 压 水 室 断 面 进行 编号 ， 
如 图 4 所 示 。 


4 叶轮 与 压 水 室 相对 位 置 示意 
Fig. 4 The relative position of the blades and the casing 
设计 流量 工 况 下 , 叶 高 中 截面 上 叶轮 区 域 的 相对 
速度 及 压 水 室 区 域 绝 对 速度 ， 如 图 5、 图 6 所 示 。 


A 


t=2T/5 


T 


201711.00626v1 


chinaXiv 


图 5 


5 设计 流量 下 叶轮 内 的 相对 速度 及 流 线 分 布 (单位 : m/s) 


Fig. 5 Distribution of relative velocity and streamlines in 


the impeller passage (unit: m/s). 


对 速度 增 大 ， 最 大 值 为 6.7m/s， 最 大 速度 位 置 发 生 在 
靠近 压 水 室 隔 舌 附近 沿 叶轮 旋转 方向 的 叶片 压力 面 出 
口 处 ; 不 同时 刻 ， 
轮 内 的 相对 速度 分 布 不 同 ， 在 局 


图 5 可 以 看 出 ， 随 着 流量 的 增加 ， 叶 轮 内 的 相 


随 着 叶片 相对 隔 舌 的 位 置 改变 ， 叶 
部 流 道内 相对 流 线 偏 


折 较 小 。 


图 6 
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6 设计 流量 下 压 水 室内 绝对 速度 及 流 线 分 布 (单位 : m/s) 


Fig. 6 Distribution of relative velocity and streamlines in 


casing passage (unit: m/s). 


勾 ， 
这 主要 是 由 于 设计 流量 工 况 下 叶片 出 口 绝对 速度 与 压 


设计 流量 工 况 下 ， 压 水 室内 的 绝对 速度 分 布 较 均 
最 大 绝对 速度 为 3.7m/s， 没 有 出 现 较 大 的 高 速 区 ， 


水 室 的 进口 绝对 速度 比较 接近 , 不 会 因 出 现 较 大 的 速度 


梯度 而 产生 冲击 ; 由 于 静止 隔 舌 的 影响 ,在 第 0 至 第 工 
断面 之 间 的 压 水 室 流 道内 出 现 低速 区 , 在 压 水 室 出 口 位 
置 出 现 较 大 的 高 速 区 。 


t=0 t=T/5 2T/5 


图 7 
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t= 4T/5 


t=3T/5 


-20.3 -17.7 -15.0 -12.3 -9.6 -6.9 -4.3 -1.6 1.1 3.8 6.4 


设计 流量 下 叶轮 及 压 水 室内 的 静 压 分 布 〈 单 位 : kPa) 


Fig. 7 Distribution of relative velocity and streamlines in 


casing passage (unit: kPa) 


压 分 布 不 同 。 整 体 上 来 看 ， 
布 相 对 均匀 ， 局 部 出 现 较 大 压力 ， 是 由 于 旋转 叶片 与 静 
止 隔 舌 的 相互 作 / 


设计 流量 工 况 下 , 不 同时 刻 叶轮 及 压 水 室 流 道内 表 
压 水 室内 及 出 口 断面 压力 分 


j 造 成 的 。 相 对 于 小 流量 工 况 ， 其 出 口 


Wt 


ACh AY ARB, (MA 1.1kPa。 


3.2 搏动 流 工 况 下 血液 汞 的 非 定常 流动 模拟 


. 能 量 计算 结果 
a 


见 图 8 。 
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8 ”搏动 流 工 况 下 泵 的 扬程 随时 间 变 化 关系 

Fig. 8 Relation between the pump head and time 

可 以 看 出 , 随 着 入 口 流量 的 周期 性 变化 , 泵 的 扬 
生 捕 动 ， 实 现 捕 动 流 工 况 ; 扬程 总 


体 变化 趋势 与 质量 流量 变化 相反 ， 流 量 较 大 的 时 刻 ， 


其 扬程 较 低 ， 扬 程 范 


典型 流量 工 况 测 点 如 图 9 所 示 。 


围 在 1.37~2.58m 之 间 。 

2. 内 流 场 模拟 结 

一 个 搏动 周期 内 泵 输出 流量 随时 间 的 变化 曲线 及 
t =0.1094s 为 整个 捕 


动 周期 的 波峰 时 刻 , 此 时 流量 接近 设计 流量 ; t, = 0.19s 
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为 整个 捕 动 周期 内 的 平均 质量 流量 时 刻 ; t, =0.2674s 
为 波 谷 时 刻 ， 此 时 质量 流量 接近 零 。 
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图 9 一 个 搏动 周期 内 的 质量 流量 变化 及 测 点 分 布 


Fig. 9 Flow rate and he monitoring points in a pulsating 


Cycle 


通过 计算 得 到 一 个 搏动 周期 内 三 个 不 同时 刻 
(t=ty ty b) 的 叶轮 区 域 的 相对 速度 分 布 ， 如 图 
10 所 示 。 


99G 


图 10 三 个 流量 时 刻 叶 轮 区 域 的 相对 速度 及 流 线 分 布 m/s) 
Fig. 10 Distribution of relative velocity at three different 


moments in impeller passage (unit: m/s). 

从 三 个 不 同时 刻 叶轮 内 的 相对 速度 分 布 可 以 看 
出 ，t 时 刻 ， 即 进口 质量 流量 接近 设计 流量 时 ,叶轮 
内 流 道 内 相对 速度 分 布 较 均匀 , 流 线 相对 光滑 t, 时 
刻 ， 流 道内 出 现 低 速 区 ,灌流 现象 明显 , 流 线 分 布 混 
AL, MATE HBL; 时 刻 , 流 道内 出 现 大 面积 
低速 区 , 每 个 流 道内 均 有 较 大 的 流 涡 , 严重 阻塞 流 道 
势必 带 来 较 大 损失 。 

图 11 为 搏动 流 工 况 下 三 个 不 同时 刻 
CGf= 丰 ty 6) 压 水 室内 的 绝对 速度 分 布 。 
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图 11 压 水 室内 的 绝对 速度 及 流 线 分 布 〈“ 单 位 : m/s) 


Fig. 11 Distribution of absolute velocity at three different 
moments in the casing (unit: m/s). 

加 时 刻 ， 即 接近 设计 流量 工 况 时 刻 ， 压 水 室内 的 
绝对 速度 分 布 均匀 ,流动 平稳 ， 流 线 光 滑 ， 此 时 流动 
损失 最 小 ; zt, 时刻， 隔 舌 附近 速度 分 布 不 均匀 ， 压 水 
室 进口 位 置 局 部 速度 较 大 , 出口 段 出 现 低 速 区 , 流 线 
也 发 生 偏 折 ; ,时刻 ， 压 水 室内 流速 分 布 不 均匀 ， 出 
口 有 较 大 的 低速 区 , 流 线 偏 折 较 大 , 并 有 较 大 的 流 涡 。 

可 以 看 出 ,搏动 流 工 况 时 ,血液 泵 的 内 部 流动 更 
加 复杂 、 均 匀 性 变 差 。 在 进行 血液 泵 优化 设计 时 ， 应 
重点 改善 小 流量 时 刻 叶轮 内 流 态 及 压 水 室 出 口 位 置 
的 流动 滞留 情况 。 

搏动 流 工 况 下 ， 三 个 不 同 流 量 时 刻 
(t=t, ty 4) 叶轮 及 压 水 室内 的 静 压 分 布 如 图 12 
所 示 。 图 13 给 出 了 三 个 时 刻 叶轮 及 压 水 室 表面 的 剪 
切 应 力 分 布 情况 ; 
叶轮 及 压 水 室 的 静 压 分 布 表明 , 不 同时 刻 , 叶轮 
及 压 水 室内 的 静 压 分 布 差别 较 大 。t 时 刻 , 叶轮 内 静 
压 从 进口 到 出 口 逐 渐 增 加 ， 压 水 室内 静 压 分 布 均匀 
亡 时 刻 ， 叶 轮 出 口 位 置 压力 分 布 不 均匀 ，, 在 压 水 室 隔 
舌 附近 出 现 局 部 低压 区 ; 4 时刻， 叶轮 及 压 水 室内 
压力 分 布 非 常 不 均 ， 将 严重 影响 泵 的 性 能 。 
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Fig. 12 The distribution of the static pressure at three 


different moments (unit: kPa). 
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13 ”叶轮 及 压 水 室 表面 的 剪 切 应 力 分 布 〈 单 位 : Pa) 
Fig. 13 Distribution of wall shear stress at three different 


moments (unit: Pa) 

图 14, 图 15 分 别 为 搏动 流 工 况 下 一 个 脉动 周期 
内 ， 叶 轮 、 盖 板 及 压 水 室 表 面 的 最 大 壁面 剪 切 应 力 
Tmax 与 平均 壁面 剪 切 应力 av 随时 间 的 变化 曲线 。 
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图 14 一 个 脉动 周期 内 叶轮 及 压 水 室 的 最 大 壁面 剪 切 应 力 


Fig. 14 Distribution of maximum wall shear stress in a 


pulsation period. 
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图 15 ”一 个 脉动 周期 内 叶轮 及 压 水 室 的 平均 壁面 剪 切 应 力 


Fig. 15 Distribution of averaged wall shear stress in a 


pulsation period 

通过 计算 得 到 t 、t, 、t 时刻 的 压 水 室内 的 最 大 壁 
面 剪 切 应 力 分 别 为 278Pa、377Pa、360Pa， 同 非 搏动 流 
工 况 类 似 , 最 大 壁面 剪 切 应 力 出 现在 隔 舌 和 叶片 进口 
出 口 位 置 。 压 水 室内 的 最 大 壁面 剪 切 应 力 高 于 叶轮 ， 而 
平均 壁面 剪 切 应 力 正 好 相反 ; 平均 壁面 剪 切 应 力 与 质量 
流量 的 变化 趋势 相同 ， 大 流量 时 刻 ， 其 平均 壁面 剪 切 应 
力 亦 较 大 ; 虽然 压 水 室内 的 最 大 剪 切 应 力 较 高 ， 但 平均 
壁面 剪 切 应 力 较 低 ， 即 整体 上 壁面 剪 切 应 力 较 小 。 在 进 
行 压 水 室 优 化 时 , 可 考虑 降低 压 水 室 局 部 的 最 大 壁面 前 
切 应 力 值 。 


4 结论 

1) 将 人 体 心脏 的 搏动 流 工 况 运 用 于 血液 汞 ， 分 析 
搏动 流动 对 无 轴 双 吸 式 血液 泵 运行 稳定 性 的 影响 , 为 血 
液 泵 在 人 工 心脏 、 人 工 透析 等 现代 机 械 辅助 系统 的 应 用 
提供 有 益 的 参考 。 

2) 搏动 流 工 况 下 ， 可 实现 泵 流量 与 扬程 的 波动 输 
出 。 一 个 脉动 周期 内 ， 各 测 点 的 压力 变化 趋势 与 泵 送 的 
质量 流量 有 关 , 相对 于 压 水 室 , 叶轮 内 的 压力 波动 较 大 。 

3) 一 个 脉动 周期 内 ， 压 水 室内 的 最 大 壁面 剪 切 应 
力 高 于 叶轮 ， 但 平均 壁面 剪 切 应 力 较 低 。 在 进行 压 水 室 
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优化 时 ， 可 考虑 降低 压 水 室 局 部 的 最 大 壁面 剪 切 应 力 
值 。 
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